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The cumulant generating function of time-averaged current is studied from an operational view- 
point. Specifically, for interacting Brownian particles under non-equilibrium conditions, we show 
that the first derivative of the cumulant generating function is equal to the expectation value of the 
current in a modified system with an extra force added, where the modified system is characterized 
by a variational principle. The formula reminds us of Einstein's fluctuation theory in equilibrium 
statistical mechanics. Furthermore, since the formula leads to the fluctuation-dissipation relation 
when the linear response regime is focused on, it is regarded as an extension of the linear response 
theory to that valid beyond the linear response regime. The formula is also related to previously 
known theories such as the Donsker-Varadhan theory, the additivity principle, and the least dissi- 
pation principle, but it is not derived from them. Examples of its application are presented for a 
driven Brownian particle on a ring subject to a periodic potential. 
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I. INTRODUCTION 

Entropy S* as a function of extensive variables charac- 
terizes macroscopic properties of materials in a unified 
manner Its functional form for a given material is 
completely determined by thermodynamic measurements 
of the heat capacity, the compressibility and the other rel- 
evant susceptibilities. Let X be an unconstrained vari- 
able of a composite isolated system. We denote the X 
dependence of S by S[X). Then, the equilibrium value 
of X is determined by 



X* = argmax[S'(X)], 

X 



(1) 



which is an example of the variational principle of ther- 
modynamics. Furthermore, according to Einstein's fluc- 
tuation theory [l|, the probability density of X is ex- 
pressed as 



Prob(X) 



^{S{X)-S{X,))/k^ 



(2) 



where fee is the Boltzmann constant, and ~ represents 
the asymptotic equivalence of the exponential factor in 
the thermodynamic limit. More precisely, let N be the 
number of the degrees of freedom. Then, the scaled en- 
tropy s{x) = limAr_j.oo S{X)/N with x = X/N fixed sat- 
isfies 



lim ^logProb(x). (3) 



In the probability theory, the right-hand side provides a 
large deviation function (up to the minus sign) That 
is, the scaled entropy s{x) is equivalent to the large de- 
viation function. 

Now, suppose that the scaled entropy s{x) is defined 
by ([3]) with the principle of the equal weight for micro- 
scopic states. Then, the variational principle ([T]) can be 
derived from the definition. A surprise in this approach 



is that the entropy thus defined characterizes thermo- 
dynamic properties completely through the fundamental 
relation of thermodynamics. After establishing the con- 
nection to thermodynamics, one can obtain the scaled en- 
tropy by thermodynamic measurements. In other words, 
the large deviation function, which characterizes macro- 
scopic fiuctuations, is obtained without measuring fluc- 
tuations. Such formulas are not restricted to the case 
of fluctuations in isolated systems. As demonstrated for 
a simple model in Appendix [XI several formulas are de- 
rived. One example is that a scaled cumulant generating 
function of a thermodynamic variable, which is given by 
the Legendre transform of the large deviation function, 
is related to a free energy density in thermodynamics. 
From this relation, it is found that the derivative of the 
scaled cumulant generating function with respect to its 
argument is equal to the expectation value of the ther- 
modynamic variable in a modifled system with an extra 
external force added. See (jA9p for such an example. The 
relations between large deviation functions and thermo- 
dynamic functions are thought to be well-known, and 
the relations lead to fluctuation-response relations, one 
of which, connects the susceptibility and the intensity of 
fluctuations of magnetization 1]. 

The main difference of non-equilibrium systems from 
equilibrium systems is the existence of a current. In par- 
ticular, the expectation value of the current in a linear 
response regime is connected with the intensity of fluctu- 
ations of time- averaged current in equilibrium 0] . This 
is one of fluctuation-dissipation relations, which might 
correspond to fluctuation-response relations in equilib- 
rium statistical mechanics. Since the time-averaged cur- 
rent appears to be important at least for small fluctua- 
tions in the linear response regime, it is natural to con- 
sider the large deviation function of time-averaged cur- 
rent. Then, from the analogy with thermodynamics, one 
may expect that there is a formula that connects the large 
deviation function with a macroscopic quantity and that 
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the fluctuation-dissipation relations are derived from the 
formula when the linear response regime is focused on. 
Indeed, Onsager and Onsager-Machlup attempted to ex- 
press the large deviation function of time-averaged cur- 
rent in terms of dissipation functions 0-0|' pointed 
out by Oono jH. This is successful for small fluctuations 
in the linear response regime, but the expression cannot 
be applied to general cases with large fluctuations. 

In the last two decades, the understanding of the large 
deviation function of time-averaged current has been sub- 
stantially developed. One remarkable result is the discov- 
ery of the fluctuation theorem which represents a simple 
and beautiful symmetry property of the large deviation 
function of time-averaged entropy production The 
symmetry property was understood from the microscopic 
reversibility of the system and it was generalized to iden- 
tities for any quantities [lol - [l^ . This progress enables us 
to revisit a linear response theory and to re-derive known 
formulas including a non-linear response formul a |15l| and 
an expression of steady state distribution [l^, Il7| . See 
Ref. [18] as a review of these re-derivations. In addition 
to the re-derivation of the known formulas, the symmetry 
property plays an essential role in deriving an extension 
of the Clausius relation (l9l - [2l| on the basis of a neat 
expression of steady state distribution [l^. As another 
direction of the progress, the large deviation function of 
time-av erag ed current was explicitly calculated for lattice 
models (23l [2^ . The results evolved into a variational 
form, which was called an additivity principle p5l - l28l |. 
When the system behavior is assumed to be described by 
fluctuating hydrodynamics, the condition for the validity 
of the additivity principle can be derived from the analy- 
sis of the path-integral expression i29h33i] . Note that the 
additivity principle is less universal than the fluctuation 
theorem, but still characterizes a universal aspect of the 
large deviation functions of time-averaged current in a 
wide range class of systems. See Ref. [34] for a recent 
related study. 

Based on these developments, we consider a possibil- 
ity that the large deviation function of time-averaged 
current is expressed by a time-averaged quantity as if 
the framework corresponds to Einstein's fluctuation the- 
ory. This formula, if it exists, is regarded as an exten- 
sion of Onsager's theory to that valid even for large fluc- 
tuations. Although the motivation seems natural, the 
problem would not be solved easily. Nevertheless, since 
a new operational variational principle associated with 
the large deviation function of time-averaged current has 
been found recently fssj , it is a good opportunity to study 
the problem. 

In this paper, specifically, we study interacting Brown- 
ian particles under non-equilibrium conditions. The mo- 
tion of the particles are assumed to be described by a 
set of Langevin equations. We then consider the large 
deviation function of time-averaged velocity of the parti- 
cles. The main results in this paper are formulas (I26p and 
(1^5)1 . The formula I^E^ indicates that the scaled cumu- 
lant generating function, which is given by the Legendre 



transform of the large deviation function, is determined 
from a variational function in a modified system with 
an inhomogeneous external force added. Essentially the 
same formula was reported in Ref. [HI . The formula ([25)1 
claims that the derivative of the scaled cumulant generat- 
ing function is equal to the time-averaged velocity in the 
modified system, which corresponds to (|A9[) in equilib- 
rium statistical mechanics. The formula (1281) also turns 
out to be an extension of the fluctuation-dissipation rela- 
tion. Since these new formulas take simple forms, we ex- 
pect that our results advance the understanding of large 
deviation function of time-averaged quantity. 

This paper is organized as follows. In section |lll we 
introduce a model and define quantities we study in this 
paper. In section Hill we first summarize main results of 
the paper and after that we derive them. In section HVl 
we address three remarks on the variational function that 
appears in our main results. Concretely, we re-express 
the variational function in terms of entropy production 
rates, mention a similarity with the Donsker-Varadhan 
formula, and explain a connection to the principle of the 
least dissipation of energy. In section |Vl we present a 
few applications of our formula by restricting our study 
to a single Brownian particle on a ring. First, on the 
basis of our formula, we re-derive exact formulas of the 
diffusion constant and the mobility in a compact manner. 
Second, starting from our formula, we show a variational 
expression associated with an additivity principle. Third, 
we propose a method for experimentally obtaining the 
scaled cumulant generating function. The final section is 
devoted to concluding remarks. Some technical details 
and related subjects are discussed in Appendices. In the 
argument below, the Boltzmann constant /cb is set to 
unity. 



II. MODEL 

We consider N colloidal particles immersed in a solvent 
at temperature T. We describe the system by a set of 
positions of the particles in a d-dimensional region f2. For 
mathematical simplicity, we impose periodic boundary 
conditions. As a typical example, one may consider a 
three-dimensional region, where one direction is confined 
by walls and no walls in other directions. Even such a 
case is included in our model by setting up a periodic 
potential whose height is extremely large at the walls. 
Another example of the region f2 is a one-dimensional 
circuit which can be generated by an optical potential. 

We denote the positions of all the particles by a; = 
{xi)f[fi, where X(fc_i)£;+; represents the Z-th component 
of the position of the fc-th particle. In this paper, bold- 
face fonts (e.g. x) represent A^d-dimensional vectors. We 
also express the position of the fc-th particle by the d- 
dimensional vector fk = {x(^k-i)d+i, X(k-i)d+2, -^Xkd)- 

The force acting on the particles consists of three types 
of force, (i) an external force, (ii) an interaction force ex- 
erted from other particles, and (iii) a force from solvent 



3 



molecules. In (i), the external force is further divided 
into a homogeneous part / = ifi)f=i and a periodic part 
determined by a potential X^fcLi Ukirk). Since the exis- 
tence of / makes the system out of equilibrium, we call / 
a driving force. In (ii) , the interaction force is determined 
by a two-body interaction potential Uk,i{rk ~ n) between 
the fc-th particle and the l-th particle under the assump- 
tion that three- and higher-body interactions among par- 
ticles are not taken into account. The total potential 
U{x) is written as 

N ^ N 

u{x) = u^^^^) + o E E - n). (4) 



fc=l 



k=l l^k 



The total force F{x) that originates from (i) and (ii) is 
expressed as 



(5) 



In (iii), the force from solvent molecules is divided into 
two parts, a friction force and a fluctuating force, where 
we assume that hydrodynamic interactions are ignored. 
The friction force is proportional to the particle veloc- 
ity Xi = dxi/dt so that it can be written in the form of 
—±iji with friction constants (7^)^'^!. We assume that 
the fluctuating force is described as Gaussian white noise 
satisfying the fluctuation-dissipation relation of the sec- 
ond kind. That is, it is expressed by ^/2T^i^i{t), where 
£,i{t) satisfies 



(6) 



(7) 



From these arguments, we obtain the equations of mo- 
tion for colloidal particles. By neglecting inertial effects, 
these are written as 



li V 7» 



(8) 



where i = 1,2,- •• , Nd. The equations are called 
Langevin equations [36| and have been studied in many 
situations including Brownian motors [37j . polymer dy- 
namics 38] and dynamics of mesoscale objects in solids. 
See Introduction of Refs. [H SO] and Chap. 11 of Ref. 
[isj . It should be noted that some recent discoveries 
of non-equilibrium relations have been tested in exper- 
iments corresponding to ([5]) with N = d = 1 [4l| - |43| . 

We assume that the particles obey a probability den- 
sity Po{x) at t = 0. After that {t > 0), the particles 
move according to (|S]). We denote by (/) the expectation 
value of / with respect to the initial distribution Po{x) 
and history of ^i{t). We define the probability density 
of x{t) = a; by P^{x,t) = {S {x{t) — x)) . Here, the su- 
perscript F of P^{x,t) represents F{x) in ([51). We use 
this superscript throughout the paper, because the force 



dependence plays an important role in our formulation. 
The evolution equation for P^{x,t) is derived as 



dP^{x,t) _ 
di ~ 

with a probability current 



Nd 

E 

i=l 



dx. 



P^ix,t) 

It 



F^{x) 



T dP^{x,t) 

li dxi 



(9) 



(10) 



The equation is called the Fokker-Planck equation [36| . 

The stationary distribution and current, which are de- 
noted by P^{x) and J^^^{x), respectively, are deter- 
mined by 

JU^)-^p^(-)--^^^ (11) 



7i 



d 



Nd 



(12) 



The expectation value of a quantity A [(a;(s))*^g] in the 
steady state is denoted by (^)^. Mathematically, the 
expectation value is calculated by replacing initial distri- 
bution Po{x) by P^{x). 

In this paper, we particularly focus on statistical prop- 
erties of time-averaged velocities 



V{t) 



dtx(t). 



(13) 



The scaled cumulant generating function G^{h) of the 
time-averaged velocities is defined by 

G^{h)^ hm ilog(e^S"-''-^-M), (14) 



where h = {hi)fj-y. In the expansion form 

r'P(h\ \ " \ " \ " "-1 "2 ■■■"'Nd 



fci=0/c2=0 fcjvd=0 



ki\k2\-.-kNd\ 



(15) 

the coefficient CE^ i. u is related to the cumulant 

ki,k2,. ..,KNd 

{Vi^'{t)V^^^{t)...V^'^-'{t))c Concretely, it can be 

shown that 

C'f,,fc.,...,fe„. = (F^(r)F2'=Hr)...y^-(r))er'=-i. 

(16) 

Here, by substituting ^ into ([T^ and taking the expec- 
tation value, we obtain 

(F,(t))=/ dx-F,{x)(- r dtP^{x,t)]. (17) 
^o" 7i \T Jo ) 



This leads to 



ml^l dx-F,(x)Pl{x). 

Jn« li 



(18) 
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Furthermore, by substituting pT|) into (IT8|) and by elim- 
inating boundary terms due to periodic boundary condi- 
tions, we obtain 



(19) 



((23l) . the expression (|23l) is equivalent to the fluctuation- 
dissipation relation (HI]). 

Now, we seek for an extension of ((23|) to that valid for 
any h and /. Here, one may recall the fluctuation the- 
orem as an extension of the fluctuation-dissipation rela- 
tion, which claims 



III. MAIN RESULTS 

The main purpose of this paper is to find an expres- 
sion of (h) which shares a common structure in equi- 
librium statistical mechanics. As reviewed in Appendix 
1X1 the cumulant generating function of a thermodynamic 
variable is related to a free energy function and this rela- 
tion leads to the important result that the first derivative 
of the cumulant generating function is equal to the expec- 
tation value of the thermodynamic variable in a modified 
system. See (|A9[) for such an example. We want to have 
a similar expression for G^{h). 

Since the relation containing (IA9p leads to the 
fluctuation-response relation for thermodynamic vari- 
ables, we expect that there might be a formula similar 
to (jA9[) . which leads to fluctuation-dissipation relations 
for currents. In order to investigate this possibility, we 
focus on the linear response regime. First, by noting 
{dG^ih)/dh,) W=o = {V.)ft, we have 



dh, ~ 



Nd 



d'^G^ih) 



0{h?). 

(20) 



dhidhj 

Here, the fluctuation-dissipation relation is written as 
'd^G^(h) 



ml 



^ Nd 

2r ^ 



dhidhj 



h=a,f=o 



(21) 

u 



See Refs. [44|, |45| for a derivation. See also Refs. 
[Tsj for a recent understanding of the derivation on the 
basis of the fluctuation theorem. The substitution of (PT|) 
into the flrst term of the right-hand side of (I^Hl) and the 
evaluation of the second term at / = yield 



dG^jh) 
dhi 



Nd 



1 

2r^ 



if,+2Th,) 



d^G^jh) 

dhidhj 



h=0,f=0 



ou\h\fh). 



(22) 



By using pT|) again in the right-hand side of p2|). we 
obtain 



dG^jh) 

dh. 



0{h\fjh). (23) 



This implies that when / and h are small, the first deriva- 
tive of G^{h) is equal to the expectation value of the 
time-averaged velocity in a modified system with an ex- 
tra force 2Th added. Since (PT|) is derived from (1^01) and 



(24) 



G^{h)^G^{-h- f/T) 



in the system under consideration [1J|. Although ([21]) 
is derived from (1^^ . a clear extension of is not ob- 
tained as far as we attempted. Putting aside the fluctu- 
ation theorem, we consider the problem from a different 
direction. The basic idea is to consider a modifled system 
with a spatially inhomogeneous extra force w{x) added. 
Precisely writing, the modifled system is described by the 
Langevin equation ([8]) with the replacement of F{x) by 
F{x) + w{x). We then define a functional ^^{w) by 



Nd 



1 



4r7, 



/ 2\P+™ 



(25) 



The first result of this paper is a variational expression 
G^(/i) =max$^(t(;). (26) 

w 

Let w^'"^^ be the optimal w that maximizes (1251) for 
given h and F. That is. 



wf^'°^^ = argmax [$^(11;)] . 



The second result is a formula 



OG^jh) 
dhi 



(27) 



(28) 



This is thought to be the generalization of ((25)) . Indeed, 
there exists a positive function iph{x) such that 



w 



d 



^ ix) = 2Th + 2T—logMx) (29) 



with a condition 



lim -— logV'h(a;) = 0. 

h-^O dXi 



(30) 



Before deriving the main results, we confirm that the 
fiuctuation-dissipation relation (j21l) is obtained from our 
resuh dH) with (1291). We start with 



dG^jh) 

dhi 



h2Th+2T 



log V'h. 



(31) 



From a basic property that there is no persistent current 
in systems without a driving force, we have an identity 



(Vi) 



log tPh. 



0. 



(32) 
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This leads to 



Nd 

~ > lini 



"5^+''-^ — 33 — 



{f,+2Th,) (33) 



-0((/ + 2T/i)2). 



By using pop . we further rewrite the right-hand side of 
(IMll as 



y Hm 



fi_dU 

d{Vi)L 
df' 



Thus, we arrive at 

Nd 



{fj+2Thj)+0{f\h\fh). 



dG^jh) 

dhi 



= V lim 



9/, 



(/, + 2Th,) 



(34) 



(35) 



0{f,h\fh). 



By differentiating ([55]) with respect to hj and setting 
f — and h. = 0, we obtain 



= 2T hni 
/=o,h=o /'^o 



dhidhj 
This is equivalent to (|2T|) . 



(36) 



A. Derivation 

The set of displacements of particles during the time 
interval lO.t] is written as 



X{t) = / dt'x{t'), (37) 
Jo 

where it should be noted that 

X{t) = V{t)t. (38) 
By differentiating (|37| with respect to t, we obtain 



^ = -F,ix) + J —i,. (39) 
at 7i y 7i 



We denote the joint probability density of x{t) and X{t) 
by p{x, X, t). The idea of considering such a joint prob- 
ability density of x and X (or accumulated work, heat) 
was employed in Re fs. |48l450l | . From the Langevin equa- 
tions ([8|) and ((39|). we derive the Fokker-Planck equation 



dt ~ ^' 



(40) 



where the Fokker-Planck operator C^p'^^ is defined by 



(x.X) 
FP -V 

Nd 



E 



1=1 



_d_F\{x) 
T 



7, ax;^" 7,9^2 

Here, we consider a quantity 

qh{x,t)^ j rfXeS-i'''^'p(a;,X,i), 
which obeys 

dqh „(x) 

-h '^h, 



7i dxidX, 



(41) 



dt 



(42) 



(43) 



where the operator Cj^ is calculated as 

Nd r / ^ -r^ , \ T 

li dxl 



E 



g F.,{x) 
dxi 



F,{x) T 2r a ■ 
hiif -I h^(p tii- — If 

li li li OXi 



(44) 



Then, from ^ and gl]), we have 

r-L 



dxqh{x,t). 



(45) 



Let /iQ be the largest eigenvalue of £[f . (143]) and (|45 
lead to the asymptotic form 



(46) 



when i is sufficiently large. Thus, from the definition of 
G-^{h) given in we obtain 



G^(M = /io- 



(47) 



The result that a scaled cumulant generating function is 
equal to the largest eigenvalue of an operator is standard 
I4I . The fluctuation theorem was studied with the aid of 
this result [13] • Recently, a relation similar to (|T7)) was 
employed for calculations of the large deviation function 
of entropy production [i^ . where a kink in the large de- 
viation function has been pointed out. See Ref. [4 71 ] for 
a recent related result. 

(x) (x)i 

We denote the adjoint operator of by , which 
is calculated as 



Nd 

E 



F^jx) d 
li dxi 

T 



T 9^ 

li dxf 



FJx), T,2 2T, d 
hiif H fij^ip-i hi- — ip 

li li li OXi 



(48) 



Let /xjj; and (n = 0,1,2,...) be all the eigenvalues of 
respectively. We label the eigenvalues 



and £^ 



6 



such that Re(AiJj) > Re(AiJ^) and Re{r]^) > Re{r]^) for 
n < m, where Re((/3) represents the real part of (p. We 
may set in general 



(49) 



with n = 0, 1, 2, where ip* represents the complex con- 
jugate of if. As shown in Appendix C, the largest eigen- 
values fj.Q and 77g are real, and thus 



(50) 



Furthermore, the eigenfunctions of the largest eigenvalue 
are positive. From (|T7)) and ([501) . we obtain 



(51) 



where iph{x) is the positive eigenfunction corresponding 
to the largest eigenvalue t^q . Here, by noting logiphix) e 
R, we define a force u^{x) as 

ul,{x) = 2T(^-^logMx) + K^ , (52) 

which corresponds to the Cole-Hopf transformation [5l| . 
By using trivial identities 



dxi 



dxi 



(53) 



2 = Vh{x) < 



we rewrite (l5Tt as 



c)l0gV'h(3;) 



F^{x)^p 



li 



uiiix) 



(54) 



(55) 



27i 7i 



Now, we consider a modified system with an arbitrary 
external force w{x) added. The stationary probability 
density in the modified system, P^^'^{x)^ satisfies PT|) 
and with the replacement of by -I- ii;. By mul- 
tiplying ((55|) by Pgf+'"(a;), integrating it with respect to 
X, and integrating by parts in the third term, we obtain 



, » Nd 



'{x) 



F.{x) 



1 



+P^^+-{^x)^uiAx) 



li 

TdP[,+-{x) 



dxi 



(56) 



Then, we substitute ((TT]) into ([55]) and we replace 
P^,+-{x)F,{x)h, - (T/TO^P.f+'^/ax, by J,^+-(a;) - 



P^{x)wi{x)/"fi. Furthermore, by substituting ([52)1 into 
the term i{x)J^~l'^ (x), we rewrite ((55)) as 



G'^(/i) = / da; V 

„ Nd jjF+w 



{x)- 



P. 



F+w 



.{xf 



(x) 



'u^,i{x) - w^{x)y 



(57) 



where we have used ([T2|) . Finally, by recalling ([191) ^-^id 
(1251) . we arrive at 



Nd 



G-(/^) = <(«.)+;^^(« 



(58) 



It should be noted that the external force w{x) in (|58|) 
is arbitrary. 

Here, the second term of the right-hand side in ([55]) is 
not negative and becomes zero only when w = u^. This 
leads to 



G^{h) = <«) > <i>^iw), 



(59) 



where the equality holds only when w — u^. This re- 
sult is the variational expression (pS)) . The optimal force 
wf^'°P\x) defined by ^ is now derived as 

w^'°^\x) = u^{x). (60) 
Next, the derivative of ([58| with respect to hi yields 



dhi 



F+w 

St 



Nd 

E 




Wi 



dh. 



F+w' 



(61) 



Since w is arbitrary, we may set w = w^'°^^ in the 
expression. This provides the result ([SS)) . Furthermore, 
from dSSl) and we obtain ^ with ([301), because the 
eigenfunction tph{x) satisfies limh^oijjhix) = const. 



IV. REMARKS ON THE VARIATIONAL 
FUNCTION $^ 

In this section, we present three remarks on the vari- 
ational function ^^{w). Before entering the main part, 
we first show a slightly different expression of ^^{w). 
Since we have (|29)) . we may restrict w{x) to the form 



d 

w(x) = 2Th- T^Vix). 

ox 



(62) 



Throughout this section, the extra force w{x) is always 
connected to the potential function V{x) through (|62p . 
We then rewrite (l26l) as 



(63) 



(/i)=max$^ (w), 
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where the optimal potential v^''°^^ satisfies 



w 



(a;) = 2Th - 



F.opt 



(64) 



Furthermore, by noting 

Nd 



Nd 



f dxJ22mj[,+-{x)= f dxY,w.ix)jf,+-{x), 

(65) 

and by using ([T^ . we express $^ (to) as 



— / dxY 



2J. 



St 



(a.) 



7« 



Wi(a;)^ 



(66) 



The expression of with was reported in Ref. 

This section is organized as follows. In section lTV A[ we 
rewrite the variational function (I66p as the difference be- 
tween two entropy production rates. In section flVBl we 
point out a similarity between the Donsker-Varadhan for- 
mula and our formula. Finally, in section flV CI we show 
that the cumulant generating function of time-averaged 
current in a simple electric circuit is expressed in terms of 
the least energy dissipation rate associated with a varia- 
tional principle that determines the voltage distribution. 



A. Expression of $^ in terms of entropy 
production rate 

In this subsection, we express the variational function 
([55)1 in terms of entropy production rate. First, the prob- 
ability current for a given probability density P{x) is 
written as 

F, , X FAx) , , T dP(x) , , 

q^{x\P) - -^Pix) (67) 

where the slightly heavy notation qf{x\P) is used in or- 
der to distinguish F and P dependence in the current. 
Here, we define the functional <j^{P) as 



By substituting (|67|) into (|68l) . we express cr^{P) as 
a^{P) 

.y^^^..g,f(.|P)(--iogP(.) + -^.(.; 



(68) 



(69) 



The integration by parts of the first term leads to 
a^iP) ^ / dxlogP{x)J2—qfix\P) 

. Nd 
+ -/ dxJ2qf{x\P)F,{x). 



(70) 



We denote the first and second terms of the right-hand 
side in ^ by i;^(P) and fl^{P), respectively. Then, 
'S^{P) is equal to the increasing rate of the Shannon 
entropy and fl^ (P) corresponds to the energy-dissipation 
rate divided by T. See Appendix [C] for a more detailed 
explanation. Since n-^{P) is identical to the increasing 
rate of the thermodynamic entropy of the heat bath, we 
identify (P) with entropy production rate of the total 
system. 

Now, from (|67p and we calculate 



F+w ^pF+u)\ _ (pF+w\ 



^-1 dxY 



St ; -^{pf^n 

Nd 

F+w („\ uF+w 



1=1 



2qf^-ix\Pf,^nM^) 



(71) 



Wi{xy 



F+w\ 



By noting q[+'^{x\P^ 
(ITTI) with (1551) . we obtain 



J^^'^{x) and comparing 



$^ (w) = i [a^+- (PJ+-) - (P^™)] . (72) 

The first term of the right-hand side in (I72p is propor- 
tional to the entropy production rate of the modified sys- 
tem in the steady state. The second term of (|72|) is a little 
bit complicated, (P^'^'^) is the entropy production 
rate just after the extra force w is turned off suddenly in 
the steady state of the modified system. In this manner, 
the expression (j72p provides us a thermodynamic inter- 
pretation of the variational function 



B. Similarity with the Donsker-Varadhan formula 

For a given path {x{t))l^Q, we define the occupation 
density as 



P™Va;) = - / dtS(x(t) - x) 

T Jo 



(73) 



We denote the probability density of {P^™{x))xen by 
Prob:^(P). We write it formally as 



Pmh^ {P) = {llS{Pr{x)-P{x)) 



(74) 



When T is sufficiently large, the asymptotic form of 
Prob^(P) is given by 



Probf (P) ~ e-^^"(^\ 



(75) 
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which is the large deviation property of the occupation 
density with the large deviation function I^{P). Let 

£pp be the adjoint Fokker-Planck operator defined by 

Nd 

^ . I .-r I n I n~ 

(76) 



E 

1=1 



FAx) d T 92 



dxi ' 7i dxf ^ 



Then, the Donsker-Varadhan formula in the Langevin 
system is written as 

(P) = - min / dxPix) / . (77) 

Such a formula was rigorously derived for general Markov 
stochastic processes [52j] . Its universal form suggests that 
it might be connected to a framework of non-equilibrium 
statistical mechanics. Indeed, in Ref. [11], it has been 
pointed out that the minimum entropy production prin- 
ciple is related to the Donsker-Varadhan formula ([77|. 
This relation was obtained by noting the detailed balance 
condition in equilibrium systems. Similarly, in equilib- 
rium dynamics, the variational principle determining the 
activity was recognized as a form similar to the Donsker- 
Varadhan formula 54]. More recently, even without the 
detailed balance condition, the variational function in the 
Donsker-Varadhan formula was conjectured to be related 
to the activity [SSj . 

Based on these achievements, we rewrite (1771) so as to 
see a similarity with our variational function (j72p . First, 
since in (|77l) is positive, we define the potential V(x) 

by 



V{x) = -2T\og(l){x). 
Then, using and we rewrite dTT]) as 



(78) 



/^(P) = — min / dxP{x) 

2T 0(x)>o Jqn 

F,{x) dVix) 



Nd 

E 

i=l 



dxi 



T d'^V{x) 



dV{x) 
dxi 



(79) 



The integration by parts in the second term leads to 

Nd 

I^{P) = --±- min / dx^ 



mm 

2T 0(x)>o Jj^jv 



F^ix) 



P[x) 



T dP{x) 
7,: 9a:,: 



i=l 

1 

27I 



dV{x) 
dxi 

dV{x) 
dxi 



P(x) 



(80) 



Thus, by using (p7|) with the replacement of F by F 
dV/dx, we obtain 



r(P) = min 



Nd 
dxY^ 



dVjx) 

dxi 



7i V 9x^ 



P{x) 
(81 



Furthermore, by comparing (jSip with (|7ip . we find 



/^(P) 



1 . 

— mm 

4 V 



dx 



F-dV/dx 



(82) 

These expressions, (|8T|) and ((82|) . were derived in Ref. 
[EE [Ell by using another method. For the large devia- 
tion function of occupation density in jump processes, a 
similar expression in terms of escape rates instead of the 
entropy production rates has been reported in Ref. (ssj . 

It is seen that there is a strong similarity between ([5^ 
and with (17^ . This suggests that two formulas 
might be described in a unified manner. We will report a 
unified treatment of the Donsker-Varadhan formula and 
ours in another paper. 



C. Relationship to the principle of the least 
dissipation of energy 

We consider an electric circuit in which n resistances 
{Ri, R2, Rn) are connected in series. We impose an 
electric potential V on the circuit by using a battery as 
shown in the left side of Fig. [TJ 



Vn 



Ri R2 



Rn 



\a] t^A 



FIG. 1: Left: n resistances {Ri, R2, Rn) are connected 
in series. We impose an electric potential V on the circuit. 
Right: By using n batteries, we impose an electric potential 
Vi on each resistance Ri. 

Let be the difference of the electric potential at the 
resistance Ri in the electric circuit. This Ki is deter- 
mined by Ohm's law 



(83) 



where is the electric current in the circuit. From the 
constraint condition — ^ ^"^^ Ohm's law 

we obtain 



V 



v„ = v 



Ri 



(84) 



(85) 



The principle of the least dissipation of energy claims 
that V<,i is characterized by the least dissipation of en- 
ergy [ssj . In order to see this statement, we consider an 
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electric circuit displayed in the right side of Fig. [T] By 
using n batteries, we impose an electric potential Vi on 
each resistance Ri such that ^ = ^- Under the 

assumption that the energy dissipation is determined by 
the Joule heating, the energy dissipation rate K is calcu- 
lated as 



y2 

K{V,,V2,...,V^) = Y^^. 

i—1 ^ 



(86) 



Then, one can directly confirm that Vn is determined as 
(KOr=i = argmini^(Fi, Fa, K), (87) 

where (X]r=i — V) means the constraint condition for 
(^)"=i- This is an example of the principle of the least 
dissipation of energy. We define the least energy dissipa- 
tion rate K^:{V) as 

K4V)= min K{Vi,V2,...,Vn). (88) 

Next, we consider fluctuations of the time- averaged 
electric current / during a time interval [0,t]. In a man- 
ner similar to (I14p . we define the scaled cumulant gener- 
ating function of the time-averaged current as 

G{h) = lim - log ieJ^^) . (89) 



For simplicity, we assume that the current fiuctuations 
are described by Johnson noise. See the first chapter of 
Ref. [36| for a brief review of Johnson noise. The Johnson 
noise satisfies the fluctuation dissipation relation 

limr((/^)-(/)^)=^. (90) 

We also assume that the fluctuations are Gaussian. Un- 
der these assumptions, we write G(K) as 



G(h) 



-h + 



i=l 1^1=1 



(91) 



Now, we connect the least energy dissipation rate 
K^iy) and the cumulant generating function G{h). 
First, from dH]) and dH]), we find 



Then, by replacing V hy V + 2Th in ([M]), we obtain 

1/2 



(92) 



K4V + 2Th) 



4TG{h), 



(93) 



where we have used (|9T|) . The subtraction of ([92|) from 
(Pg)) leads to 

G{h)^^[K,{V + 2Th) -K,{V)]. (94) 

This result that G{h) is written by the difference of the 
energy dissipation rates for two different states reminds 
us of (|A8[) in equilibrium statistical mechanics. Our for- 
mula (|63p with ((7^ is a generalization of ([M)) . 



V. APPLICATIONS 

In this section, we present three applications of our 
result by restricting our investigation to a Brownian par- 
ticle on a ring. We write an explicit setting. Let x{t) be 
a position of the particle along the ring of size L. We 
assume that the motion of the particle is described by 



7 V 7 



(95) 



with the periodic boundary condition, where F{x) is a 
deterministic force given by 



F{x) = / - 



dU{x) 
dx 



(96) 



Such a force can be realized by using a rotating periodic 
potential [Mlii]. Note that (|95]) with ^ corresponds 
to the case that TV = d = 1 in the model studied in 
Section ini The stationary distribution P^{x) and the 
stationary current are determined by 



1 



j;, = fF(x)p,^(x) 



T d p 

jdx 



(^) 



(97) 



with the periodic boundary condition P^{0) — P^{L) 
and the normalization condition dxP^{x) = 1. 

For this system, in section fVAl we derive the diffusion 
constant D = limT-_j.oo(l/2)r (^(t)^)^ and the mobil- 
ity PL = (d/df) limi-^oo {V{t)) c by utilizing our formula. 
Then, in section IV Bl we show that our variational ex- 
pression (1^51) with leads to a formula derived from 
an additivity principle [25ii - t28j . Finally, in section IV CI 
we present an experimental method for determining the 
cumulant generating function without the information of 
F{x). 



A. Diffusion constant and mobility 

In the system under study, the stationary distribution 
P^{x), the diffusion constant D and the mobility p are 
known to be expressed in terms of the following functions 

[lililisi: 



Jo 



Concretely, the expressions are 

P^(a;)=a/_(x), 



and 



T 

D ^ -L^a^ / dxl+{x)l-{xf 
7 Jo 

L 

fi = —a^ / dxl+{x)l-{x), 
7 Jo 



(98) 



(99) 



(100) 



(101) 
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where a is the normaUzation constant given by 

dxl±{x) 
/ 



(102) 



Here, one may wonder about the origin of such compact 
expressions. Below, within our formulation, we uncover 
a simple structure behind the expressions. 

First, we consider the diffusion constant. We start with 
the formula ([58]) with (l25|). In the case that N = d ^ I, 
it becomes 



h{v):r 



(-^)r-(K--)'> 



2\ F+w' 

1 

St 

(103) 



where it should be noted that w{x) is arbitrary. By set- 
ting w = in (|103p . we obtain 



G^ih)^h{V)f, + ^{{u()')^^. (104) 



Since uf{x) 
panded as 



h=0 



(See (Uni) and uf^{x) is ex- 



(105) 



n=l 



The substitution of (fTOS]) into ([T04| yields 



G^{h) ^ h {V)f, + — ((^.l)^); + 0{h'), (106) 



4T7 



which leads to 



D 



4T7 



(107) 



Next, we consider the mobility. Let us fix a function 
w'{x). By setting w — —F + w' in (|103p . we rewrite (|103p 
as 

G^{h) = h {V):y^ {{w' - Ff - « -w' + F)X ■ 

^ (108) 
We then differentiate (|108p with respect to force /. The 
result is 



dG^jh) 
df 



(109) 

Since w'{x) is arbitrary in this expression, we may set 
w' = + F. This leads to 



dG^jh) 
df 



1 / F\P'+< 



(110) 



By substituting (IT05l) . (IT06l) into plO)) . we obtain 

1 



2r7 



(111) 



As shown in (|107p and (lllip , D and are expressed by 
the second and the first moments of ui{x), respectively. 
Here, we determine ui{x). We substitute (I105|) into ([55)) 
and extract the terms proportional to h. We then obtain 

2T {V)f, = -^u,{x) + -F{x)u^{x). (112) 
7 ax 7 

Note the periodic boundary condition iti(O) — ui{L) and 
the normalization condition 



/ dxui{x) ~ 2TL, 
Jq 



(113) 



which comes from (|52l) . By comparing (|112p with (|97l) . 
we find that ui{x) can be solved in the manner similar 
to P^^ix). The result is 



ui{x) — 2TLaIj^(x). 



(114) 



By substituting ([TO]) and ^ into ([TUT)) and (ITTT)) . we 
obtain ()100p and ([TUT]). Furthermore, we can systemat- 
ically evaluate higher-order cumulant coefficients by for- 
mulating an iterative relation. See Appendix |D] 

At the end of this subsection, we study the violation 
of the fluctuation-dissipation relation. First, from ()18l) . 
we have 

{v)i^< - ^ {F)::< = {ui)T'^ . (115) 



We substitute ()115l) into (|110p and combine it with our 
main result (1^51). We then obtain 



dG^jh) _ 1 dG^jh) 
df ~ 2T dh 

which leads to 



2T-1 ^ 



^-^.^limA(f)-+< 



(116) 



(117) 



Furthermore, by considering an expansion form 

Pj+""(x) = P[,{x) + hQ,{x) + 0{h'), (118) 
we rewrite the right-hand side of (|117|) as 



1 



2^7 Jo 



dxF{x)Qi{x). 



(119) 



The expression ()119p represents a degree of the violation 
of the fluctuation-dissipation relation. Recently, several 
attempts have been presented so as to characterize the 
violation of fluctuation-dissipation relation [60-64j. It 
might be interesting to find a connection of ()119p with 
them. 
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When the driving force / is absent, Qi{x) becomes 
zero so that the fluctuation-dissipation relation holds. 
We show this fact explicitly. We substitute ()118|) into 



(P?]) with the replacement of F{x) with F{a 



4. (a;), 



and extract the terms proportional to h. We then obtain 
the following differential equation for Qi{x): 



d_ 

dx 



-F{x)Q^{x) + -u^{x)P[,{x) - --^Qiix) 
7 7 'J ax 



(120) 

with the periodic boundary condition Qi(0) — Qi{L) and 
the normalization condition 



dxQi{x) — 0, 



(121) 



which comes from (|118|) . Here, from (j99| and (|114|) . 
we find that ui{x)P^{x) is constant in x when F{x) = 
—dU{x)/dx. In this case, (|120l) provides us Qi{x) = 
Cie~^(^)/^, where d is a constant. Then, (IT2T1) leads 
to Ci =0. Thus, recalling (|119p . we conclude that the 
right-hand side of (|117p is equal to zero when / is absent. 
Note that the discussion above applies to many-body sys- 
tems and higher-dimensional cases. 



B. Relation to an additivity principle 



We start with the variational expression ((26)) with (|25|) . 
Since the probability current is independent of x in the 
one-dimensional case, we can rewrite (I26|) as 



G^(/.) 



max 



F+w 

St 



L 



dx 



F 



F+w 



St 



(^) 



4T7 



w{x)'^ 



(122) 

where J^'^^ and P^^'^{x) are the stationary probabil- 
ity current and density in the modified system. Note 
that J^'^^ and P^^'^{x) can be determined uniquely 
from (|97p for a given extra force w{x). Conversely, when 
a constant current J and a probability density function 
P{x) are given, we can define the extra force w{x) by 



w{x) 



Fix) 



TdF{x) P{x) 
^ dx 7 



F{x) 



(123) 



Thus, there is a one-to-one correspondence between w{x) 
and ( Jg^^™, Pg^^"'(a;)) . From this fact, we can rewrite 
(IT^ as 



(h) — max 
p,j 



hJL 



— — / dx 



1 



AT 



F{x) 



J 



T dPjx) _ F{x) 
J dx 7 



F{x) 



(124) 



By taking the maximum of P{x) first, we obtain 



G (h) — max 
J 



hJL 



—— mm 

4T p 



'^pM 



T dF{x) F{x) 
J dx 7 



F{x) 



(125) 



Therefore, we find that the Legendre transform of G^{h) 
with respect to hL, which is the large deviation function 
/(J), becomes 



/(J) 



— - mm 

4r p 



f 

Jo 



dx 



P{x) 



T dF{x) P{x) 
J dx 7 



F{x) 



(126) 



This expression is the same as the variational expression 
derived from the additivity principle p5l - l28| . 

When the system is described by fluctuating hydro- 
dynamics, the variational expression is obtained under 
the assumption that the most dominant contribution to 
a path integration is given as a stationary configuration 
|29l - l33{ . Here, since the Langevin equation ^ can be 
mapped to a stochastic partial differential equation [63] , 
which is equivalent to fiuctuating hydrodynamics, we can 
derive the variational expression in our system if the as- 
sumption is satisfied. Our result (|126p indicates that the 
assumption is valid. 

In the derivation of (|126p . the spatial homogeneity of 
the stationary current plays an important role. However, 
since this property is specific to one-dimensional systems, 
it is not obvious whether the same variational principle 
does hold in higher-dimensional cases |66i | . Furthermore, 
by noticing that the probability current and the probabil- 
ity density considered here are equivalent to the particle 
current and the particle density in non-interacting parti- 
cle systems, one might find that it is a highly non-trivial 
problem to derive (jl26p for interacting particle systems. 



C. Experimental determination of cumulant 
generating function 

We consider an experimental method for determining 
G^{h) without the information of F{x). Since the direct 
evaluation of cumulants becomes harder and harder as 
the order of the cumulants increases, it may be useful if 
we could present an efficient method for evaluating G^{h) 
on the basis of our formula under the assumption that we 
measure only the trajectories (x(i))^^Q and values of T 
and 7. 

The basic idea is to utilize (|58l) . which was a corner- 
stone in the derivation of the main results ([25)) and ()28l) . 
First, we add some external forces w{x) to the system 
and measure trajectories of the particle. We then notice 
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that the left-hand side of ([58)1 is independent of external 
forces w{x). Since we can estimate the right-hand side 
of (|58p by using the measured data, we obtain equations 
which yield it;^'°^'(a;). Essentially the same idea was re- 
ported in Ref. [35| • Below, we present a slightly different 
one from that in Ref. [35j |. 

Concretely, in terms of a complete set of orthogonal 
functions {4>a)^=o, we expand as 



F,opt 



ix) 



a=0 



M^(/'a(a:), 



(127) 



where m is a truncation number of the expansion. As 
m is increased, the accuracy of the expansion (jl27p in- 
creases. For later convenience, we choose 4>oix) = 1. 
Since dxw^'°'^\x) = 2ThL (See the coefhcient 

Uq is derived as 



= 2Th. 



(128) 



We shall determine the other expansion coefficients 



First, by differentiating 
tain 



with respect to h, we ob- 



1 




F,opt 



dh ^' ■ 2T7 

We then expand w{x) as 

m 

W{x) = ^ Wa(l)a{x). 



dh 



F+w 



St 

(129) 



(130) 



a=0 



By substituting (IT27|) and (fT30)) into we rewrite 

(fng)l as 



dh 



Here, by defining 



' a=Q 



{u'^ - Wa){(f>a(t)ri 



,F+w 
'st J 



~dh' 

(131) 



(132) 



we rewrite (I13ip as 
dG^ih) 



dh 



= {V) 



F+w 

St 



2TAl[ 



-jr- (133) 



where we have used (|128p . We prepare m + 1 external 
forces w''°\x),w'^'^^x),...,w'-"'^x). Recalling that 
is valid for any external force w{x), we equate the right- 
hand side of (|133p for the case w^^^ (x) to the right-hand 



side of (|133|) for the other cases ^^"^(2;) (a = l,...,m) 
The obtained equations are 



F+w 

St 



(0) (0) 

- 2TAf - {V) 



St 



2TA^ 



a: 



„(o)\ 



dh 



(134) 



where a = 1, m. We further define an m x m matrix 
as 



4 = 4"^ 



Al 



. (135) 
With the assumption that A is invcrtible, we derive 



du 



dh 



^1 = Y.{A'\^ \{V) 



1st 



2TA] 



,(0) 



(136) 



mi 



2TAI 



This expression implies that du^/dh is expressed in 



Now, 



terms of (u^J)™ 1, (y)f/'" and {^s) 
from trajectories of the particle in the modified systems 
with the external forces w^"^ (x) added, we approximately 

obtain {V) 



F+w'-' 

St 



and 



1 

T JO 



in the form 
dtx{t), (137) 



\F+w° 
/st 



dtq}„ {x{t)) (bs (xit)) , (138) 



with < rj, 6 < m. Note that the approximation 
becomes exact in the limit t — > cxd. Therefore, we 
can express du^/dh by {u^)'^^i from experimental data. 
By integrating du^/dh with the condition u^^\h=Q = 0, 
{I < rj < m), we obtain (u^)^i as a function of h. 
Similarly, (11331) leads to the expression of dG^ (h) j dh in 
terms of {u'f^)'^^i- The integration yields G^{h), where 
we have used G^{h)\ii=o — 0. 

We performed a numerical experiment of the driven 
Brownian particle on a ring subject to a periodic poten- 
tial U{x) — Uocos{2nx/L). We set ipkix) = cos{2TTkx/L) 
(0 < fc < m/2), 0fc(a;) = sin (27ra;(fc - m/2)/L) (m/2 < 
k < m). The m + \ external forces (wfe(a::))^Q are cho- 
sen as {C4sk{x))™^Q, where C is a constant. Following 
the recipe we described above, we obtain G^{h) and 
K;^'°^'(a;) from trajectories of the particle. The results 
are displayed in Fig. O We also calculate G^{h) and 
it;^'°^'(a;) by approximately evaluating the largest eigen- 
value and the corresponding cigenfunction of the opera- 
tor See Appendix lEl for the calculation method. 
These two results are close to each other. It turns out 
that our experimental determination method is useful to 
obtain the accurate value of G^(h) from measured tra- 
jectories. 
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FIG. 2: (Color online) Numerical experiment for measure- 
ment of G^{h) (up) and Wf^'°^^{x) with h = 3 (down). 
Quantities are converted to dimensionless forms by setting 
7 = r = L = 1. We fix / = 1 and C/o = 3. We assumed to 
know T = 7 = 1 and experimentally determined (h) and 



F,opt 



(x) from trajectories (a;(t))J^o following the method de- 
scribed in the text. We set m = 10, C = 10 and r = 400000. 
By taking 10 samples, we estimated (h) and w^'°'^^{x). 
The obtained results are displayed with green dashed lines. 
Error-bars are within the lines. The red lines were obtained 
from the evaluation of the largest eigenvalue of 



VI. CONCLUDING REMARKS 

In this paper, we have presented the new formula 
for the cumulant generating function of time-averaged 
velocity. The formula enables us to determine the cu- 
mulant generating function in terms of the expectation 
value of time-averaged current in a modified system. This 
corresponds to (jA9[) in the example of equilibrium statis- 
tical mechanics, and it is also regarded as an extension 
of the fluctuation-dissipation relation. Furthermore, we 
have pointed out that our formula is related to the addi- 
tivity principle, the Donsker-Varadhan formula, and the 
least dissipation principle. Before ending the paper, we 
address problems that should be studied in future. 

First, since our formulas are rather formal, we should 
explore further applications of our formulation to vari- 
ous physical problems. With regard to this aspect, in 



Appendix [Fl we present a simple model in which we can 
explicitly calculate the cumulant generating function and 
the optimal force. We hope that this model might help 
readers to grasp our formulation and promote the appli- 
cations. Furthermore, although our formulation is valid 
for interacting particle systems, useful predictions are not 
presented yet. A key step in the analysis of many-body 
problems is to extract slowly varying effective variables 
from many-body distribution functions. Here, we may 
employ some approximations for extracting effective vari- 
ables. The result of the approximate calculation with our 
formulation might provide a new insight into properties 
of current fluctuations because our formulation is quali- 
tatively new. 

Next, we ask the range of the applicability of our for- 
mulation. The first question may be to consider similar 
formulas in under-damped Langevin systems. The for- 
mulas that correspond to ([^5]) and (I^S)) can be derived 
in a straight-forward manner. This fact is not surpris- 
ing at all when the formulation of this paper is properly 
understood. However, it may be difficult to find a physi- 
cal interpretation of the formulas. For example, it seems 
that "the extra force" appearing in the variational for- 
mula cannot be realized in laboratory experiments. More 
importantly, as far as we attempted, the variational func- 
tion may not be expressed in terms of the entropy pro- 
duction rates, which is different from the over-damped 
case. We do not understand whether or not these are 
fatal in our formulation. We need further investigations. 

Here, with regard to the range of the applicability of 
our formulas, we restrict our concerns to a formal aspect 
by putting aside physical interpretations. Concretely, in 
Markov jump processes, we can derive the formulas that 
correspond to ^ and As seen in (|G14I) and (jGlSP 
in Appendix [Gl the mathematical expressions are quite 
simple. We thus expect that there might be physical ex- 
amples for which these formulas are useful. Furthermore, 
since these are expressed in terms of escape rates, it is 
conjectured that similar variational principles are formu- 
lated in many-body Hamiltonian systems. It is stimu- 
lating to prove this conjecture. By combining physical 
interpretations and formal aspects, we wish to seek for a 
universal structure behind our formulation. 

Finally, we go back to our basic motivation. Our naive 
idea is that the large deviation theory for current fluctua- 
tions might be associated with an operational framework. 
This was conjectured from the well-established fact that 
the large deviation theory for fluctuations of thermody- 
namic variables is tightly connected to equilibrium ther- 
modynamics. The progress developed in this paper does 
not provide a final answer to the problem, but the results 
are encouraging us to continue our efforts. In the best 
case, we might have an operational framework of thermo- 
dynamics for currents. Then, we will consider a possibil- 
ity to unify such a framework with equilibrium thermo- 
dynamics. It may be plausible that a unified scheme is 
related to the so-called steady state thermodynamics [l9l - 
[2ll . [67l - [69| . With dreaming such possibilities, we chal- 
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lenge to study operational viewpoints in non-equilibrium 
statistical mechanics. 
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Appendix A: Equilibrium statistical mechanics 

We review the relation between the cumulant generat- 
ing function of a thermodynamic variable and a thermo- 
dynamic function. Specifically, we consider a spin model 
with finite range interaction in a finite dimensional lat- 
tice A. Let (Ti, i € A, be a spin variable whose value is 
either -|-1 or —1. We denote {<^i)i^\ by a collectively. 
We assume that the Hamiltonian for the system under a 
magnetic field H is given by 



(Al) 



where is the magnetic moment. According to equilib- 
rium statistical mechanics, for any observable A{a), the 
expectation value {A)t.h in the system at temperature 
T is given by 



3-Wff(cT)/T 

Z{T,H) ' 



where 



(A2) 



(A3) 



As an observable, we consider the magnetization density 
'"(^) = ^E^- (A4) 



where N — |A|. Now, we define the free energy density 
f{T,H) by 



/(T,i/) = -r Jim -logZiT,H). (A5) 

Then, the fundamental relation of thermodynamics is 
written as 



df 



-sdT — mdH, 



(A6) 



where s is the entropy density and m is the expecta- 
tion value of the magnetization density m{a). Notice that 



/(T, H) is determined from the measurement of the heat 
capacity c(r, H) and m{T, H). 

Putting aside the thermodynamic function, we con- 
sider the scaled cumulant generating function of m((T), 
which is defined as 

GrAh)^ lim ^\og{e''^^''^^)T,H ■ (A7) 



From (jASP and (jA7[) . we find the following relation be- 
tween the cumulant generating function and the free en- 
ergy: 

GtmW = -^[f{T,H + hT) - f{T, H)] . (A8) 

In particular, from the thermodynamic relation (jA6p . we 
obtain 



dh 



miT,H + hT) 



(A9) 



for any h. By differentiating (|A9[) with respect to h and 
set h ~ 0. we obtain 



lim N((m(a)-mf) 

N^oo \ /t,H 



= T 



dm{T,H) 



, (AlO) 



N^oo \ ' ' ' ' / T.H dH 

which is referred to as a fluctuation-response relation [l| . 

Appendix B: Properties of the largest eigenvalue 

We show that the largest eigenvalue of the operators 
and is real and that the corresponding eigen- 

functions are positive. The argument presented here is 
not mathematically rigorous, but provides a practical 
idea for understanding the claim. Although we focus 
on one-dimensional case for simplicity, the extension to 
higher-dimensional cases are straight-forward. We ex- 
press the operator in the form 



C^f^'^ .^ = A(x) — ^ + B(x) — ^ + C{x)^, (Bl) 



where j4(a;) > 0. We study the eigenvalue problem of /^Ij^-* 
in the space of periodic functions defined on the interval 
[0,L]. 

We define Xi = Li/N, i = 1,2, N, and let Aa; be the 
distance between x^+i — Xi. We replace the differential 
operator d/dx and d^/dx^ by differences. Concretely, 
the following replacements are considered: 



dip 
dx 



9V 
dx'^ 



2Ax 



(pjxj+i) + (pjxj^i) - 2ip{xi) 
(Ax)2 



(B2) 



(B3) 
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for i = 2,3, ...,N — 1. A special care is needed for the 
cases i — 1 and i ~ N. For example, 



dip 
dx 



ip{x2) - (pjxN) 

2Ax 



5V 
dx"^ 



(fixi) + ip{xN-l) - 'ilfixN) 



X—XN (Ax)"^ 

Then, we define an TV x TV matrix M^v such that 

(MNip)i 



(B4) 



(B5) 



= A{xi) 



dx"^ 



+ B{X^) 



dip 
dx 



+ C{xi)ip{xi), 



(B6) 



where the iV dimensional vector is defined by — 
ip{xi). Since Mjv is an approximation of the operator 
C''^\ we assume that the eigenvalues and eigenfunctions 

(x\ 

of C]^ are determined from the limit iV — > cx) for the 
eigenvalues and eigenvectors of Mm- Thus, below, we 
study the eigenvalue problem of instead of the eigen- 

{x) 

value problem oi C\ . 

We show that the largest eigenvalue of is real and 
that all components of the corresponding left and right 
eigenvectors are positive, when N is sufhciently large. 
This claim is understood from the Perron- Frobenius the- 
ory. A real matrix B is called an irreducible ML-matrix 
[7(J | if B satisfies the following two conditions: (i) All 
non-diagonal components of B are non-negative, (ii) 
There exists a positive number a such that the matrix 
T = al + B is irreducible, where / is the unit matrix. 
The Perron- Frobenius theory leads to the statement that 
the irreducible ML-matrix B has an eigenvalue fipt that 
satisfies the following conditions: 

1. /ipf is real. 

2. The associated left and right eigenvectors are 
strictly positive and unique up to constant mul- 
tiples. 



probability density is given by P. We define a probability 
density p^{x, t) as the solution of 



-p^ic,,t)=-J2—qrix\p^ 

i—1 

with the initial condition 



(CI) 



(C2) 



Here, we take the limit t — > in (IC1[) . By noting (jC2l) in 
the right-hand side of (|Cip . we obtain 



f) Nd ^ 

lim-p-(.,0 = -5:-,f(.|P). (C3) 



dt ' ^ dx 

1=1 

Thus, the first term of ([7D|) . (P) is expressed as 



(C4) 



Since J^jv dxp^{x,t) — 1, (|C4[) is further rewritten as 
d 



(P) = 



dt 



/ dxp^ {x,t) log {x,t) 
Jo" 



"(C5) 

This shows that Y,^{P) is the increasing rate of the Shan- 
non entropy for a given distribution P. 

Next, we consider the second term, il^{P). We shall 
show that n^{P) is the expectation value of the energy 
dissipation rate divided by T when the distribution is 
given by P. According to the energetic interpretation of 
Langevin equations [Til, [tI] , the energy dissipation rate 
^! for each realization is given by 



Nd 



(C6) 



where the symbol o represents the multiplication in the 
sense of the Stratonovich [sB]- It should be noted that 
is a stochastic variable. By substituting (jS]) into (|C6|) . 
we write 



3. For any eigenvalue ^ except for /ipf, the real part ^ ~ 5^ ( ^^'^^^'^^ 

of /I is less than /Ltpf. 




Fi{x) o 



(C7) 



See Theorem 2.6 of Ref. [70|- Since A{x) > 0, the matrix 
Mjv for sufficiently large iV turns out to be an irreducible 
ML-matrix. We thus obtain the claim that the largest 
eigenvalue of Mm is real and that all components of the 
corresponding left and right eigenvectors are positive. 



Appendix C: {P) and entropy production rate 

We explain thermodynamic interpretations of the first 
and second terms of (|70p . which were denoted by S^(P) 
and Vl^ {P) in Section HV Al respectively. 

We first consider the first term, T,^{P). We assume 
that the system is set under the force F and that the 



Here, we change the multiplication rule from the 
Stratonovich type (o) to the Ito type (•) [11]. Then, (fCTj) 
becomes 

(C8) 

By taking the expectation value of (jC8[) and integrating 
by parts the second term, we obtain 



Nd , 
= / dxY,F^{x) -F,{x)P{x) 



T dP{x) 

li dxi 



TVt^{P), 



(C9) 
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where we have used {Fiix) ■ ^j) = 0. The last expression 
claims that Vl^ (P) is the expectation value of the energy 
dissipation rate divided by T when the distribution is 
given by P. 



Appendix D: Iterative expression of cumulant 
coefficients 



For the system studied in Section |Vl we determine the 
cumulant coefficients in an iterative manner. We start 
with the expression 



2TG'' (h) 



J dx 7 



Fix) 



which is obtained from (|55|) with d — N — 1, where 
uj^{x) satisfies 

[ dxul [x) = 2ThL, (D2) 
Jo 



ul[Q)^ul{L). 
From a trivial relation 



dx 



T 



= e- lo dyF(y)/T 

we rewrite (|Dip as 
d 



d_ 

dx 



ul{x)e^Sdyny)/T 



(D3) 



(D4) 



dx 



(D5) 



u'[{x)e^o''y^'^y^/^ 

By integrating (jD5[) with respect to we obtain 

(D6) 

where C is an integration constant. Here, by direct cal- 
culation, we can derive an identity 



dyg{y)e^o '^^n^)/^ - dyg{x + y)e^S^' "^^^^^^ 
Jo 

(D7) 



for any periodic function g{x). We then set 
g{x)=2^G^{h)-l^ul{xf 



(D8) 



in (jD7l) and substitute it into the second term of the 
right-hand side of (|D6|) . The result is 



C- 



1 

1 - ef^/^ 



u'[{x)^e-SSdvF(y)/T 

X dyeJr^ '^^^W/^ (2^G^ih) - A^.f (. + yf 

(D9) 

where C is a constant. The condition u^(0) = uf^{L) 
leads to C = 0. Through the transformation of the inte- 
gration variable y L — y, we rewrite (jD9[) as 



(Dl) u^ix) = 



27 



1 - e-f^/T 



G^{h)I+{x) 



1 



4^7 7o 



dyefr'dzFi.)/T^F(^^__yy 



(DIO) 



where we have used (|98p . By integrating the both sides 
of (IDIOI) with respect to x and using the condition (|D2I) . 
we derive 

G^{h) = hL-{l - e~f^l^)a + dxl.{x)ul{x)\ 
7 4r7 Jo 

(Dll) 

where a was defined in (|102p . 

It should be noted that (x) satisfies a self-consistent 
equation. Indeed, by substituting ()Dlip into (|D10|) . we 
find 



ui(x) =2ThLaI+{x) 



J:-ydzF{z)/T 



2T{1 - e'f^/T) Jo 



dy 



a c?r/_(r)ii^ (r) ~u^{x-y) 



(D12) 



By substituting dTUSj) into (lDT2t . we can determine the 
coefficient u„(x) itcratively. Concretely, by noting 



oo oo 



oo /n— 1 ^ 

ri=2 \i=l / 

we derive a recursion equation 

n-l „i 



(D13) 



Unix) = 



1 



2T(l-e-/i/^)^ 



/;-«d^F(z)/T 



a I / drl^{r)ui{r)un^i{r)\ - ui{x - y)Un-i{x - y) 



(D14) 
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for n > 2. For example, by recalling (|114l) . we calculate 



U2{x) 



aij^ drl_{r)l+{rf 



• dzF{z)/T 



(D15) 



Furthermore, we expand G^{h) as 

GP{h) = y ^/i". 

71 = 1 

By substituting (|DT3|) and ((DT6| into (|DTT|) . we obtain 

Gi= L-il~e-f^/^)a (D17) 
7 

and 



an! 



4^7 Jo 



Gn = J2i^. I dxI^{x)ui{x)Un-l{x) (D18) 



for n > 2. With this formula, the cumulant coefficients 
Gn can be calculated in an iterative manner. For exam- 
ple, the third-order cumulant coefficient is derived as 



G. = 



7 



dxl^ (x)I+{x)u2{x) . 



(D19) 



Reimann et al found the exact expression of the diffu- 
sion constant (|100l) from the first and second cumulants 
of first passage time (3§|. Here, we remember an iter- 
ative formula for the moments of the first passage time 
[36| . By comparing this iterative formula with our iter- 
ative formula, one might find a general relationship be- 
tween the cumulant generating function and statistical 
properties of the first passage time beyond the first two 
cumulants. 



Appendix E: Calculation of G{h) and w'^^{x) 

We focus on the case that N = 1 and d = 1, which 
is studied in section |Vl We numerically solve the largest 
eigenvalue problem (|5T|) . First, we consider a linear equa- 
tion 



d 



f^cl){x,t) ^ c[^^U{x,t). 



(El) 



Since G^{h) is the largest eigenvalue of cl^^\ the long 
time behaviors of solutions are described as 



(E2) 



Now, we assume (/>(a;, t) > for any t by choosing appro- 
priate initial conditions. Then, we define 

u{x,t) = 2T—log(j){x,t) +2Th. (E3) 



By dividing (lEll) by (f>{x, t) and using the same method 
deriving (|55|) in section III.A, we obtain 



d F(x) 1 T d 

2r— log0(a:,i) = -^uix,t)+—uix,t)^+-—u{x,t). 
ot 7 27 7 ax 

(E4) 

Finally, by differentiating (jE4|) with respect to x, we 
rewrite (IE4I) as 



(E5) 



(D16) with a current 



F{x) 



1 T d 

ju{x,t) = -^-^^u{x,t) - —u{x,tf —u{x,t). 

7 27 7 ox 

(E6) 

Let u»(a;) be the stationary solution given as the long- 
time limit of u(a;, t). Then, from (|E2I) and (|E3p . we obtain 

d 

u^x) ^2Th—\ogi!h{x) + 2Th. (E7) 
ox 

By recalling ([29|) . we have 

u.°P'(x) (E8) 
Furthermore, by comparing (jE6|) with (j55p . we obtain 

1 



G'{h)^-—\imUx,t). 



(E9) 



Therefore, by numerically solving (jESp . we can calculate 
G^{h) and wl'^ix). 

Appendix F: Exactly calculable model 

We illustrate our theoretical idea by explicitly deriving 
the cumulant generating function and the optimal force 
for a linear Langevin equation. The model describes the 
motion of a single Brownian particle in two dimensions, 
which corresponds to the case that N = 1 and d = 2 in 
(|H]). By setting 7 = 1 for notational simplicity, we write 
the equation as 



dx 



= F{x) + v^l. 



with 



F{x) = -kx + el{x), 



(Fl) 



(F2) 



where I = {—X2,xi). From (fTTjl and (fT2|) . we obtain the 
stationary distribution as 



P,^{x)^Cexp(-—x'], 



(F3) 



with a normalization constant C. Although the station- 
ary distribution is Gaussian, the detailed balance condi- 
tion is not satisfied because of the circulation force el. 
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One can check this fact immediately from (fill) and (jF3l) . 
It should be noted that the linear Langevin equation was 
studied with a novel idea, irreversible circulation, which 



characterizes the violation of detailed balance [73| . 

The quantity we consider here is a time-averaged an- 
gular momentum (or irreversible circulation) defined as 



L{r) 



dtl{x{t)) ■ x{t). 



For the scaled cumulant generating function of L, 
G{h) = lim - log /e^''-^(^) 



(F4) 



(F5) 



we can derive the variational principle by using the same 
method as Section UlIAl First, we define an operator Zh 
as 

/:„-^ = TvV-v.(F^) + /i(rF)^ 

+ h^Tfip^2ThV -{lip). 

Since G{h) is equal to the largest eigenvalue of the oper- 
ators Ch and we have 



£), • = G{h)^, 



(F7) 



where is the positive eigenfunction corresponding to 
the largest eigenvalue of C]^. Next, we define a force u as 



u{x) = 2T (V log (j3{x) + hl{x)) . 



(F8) 



By dividing (jF7p by (j) and using the definition of it, we 
obtain 

F -u + ^^v? + TV ■u = 2TG{h). (F9) 



This takes the same form as ((55|) so that we can follow 
the argument from ((55|) to ([58|) . We then arrive at 



G{h) = ^,,{w) + ^{{w-uffJ'" 



with 

^h{w) dxl • 

Jn 

Here, it should be noted that 



n 4T\ 



Jf2 



(FIO) 



(Fll) 



(F12) 



(jFlOp is the key identity of our formulation, which cor- 
responds to ((55)) . Indeed, from (jF10[) we obtain 



G(/i) = max$/i(t(;) 



(F13) 



(F14) 



where the optimal force i^opt is given as 
Wopt{x) = u{x). 



(F15) 



Finally, by solving (|F9p . we derive the explicit expres- 
sion of u and G{h). Suppose that 

= exp(-^a;2) , (F16) 

where a is a constant. This yields 

u{x) = -2Tax + 2Thl{x). (F17) 

By substituting (IF17|) into dEH), we obtain 

G{h) = -2Ta + x^ [Ta^ + ka + h{Th + e)] . (FIB) 

Since the left-hand side is independent of x, we determine 
the constant a as 



(F19) 



Here, we can confirm that (IF7p holds for (|F16p with 
(IF19I) . because ((F7| is equivalent to (|F9t . We then obtain 



G(/i) = /c - y/k^ - iThiTh + e). (F20) 



From (jF17[) and (jF19p . the force in the optimal system 
is explicitly written as 



F{x) + u{x) = -y^k^ -4Th{Th + ()x + (e + 2Th)l{x). 

(F21) 



Appendix G: Formula in Markov jump process 

The model we consider here is a continuous-time 
Markov process of a finite space M.. We assume that 
the transition rate, w(a a') for a 7^ a' e A^, is irre- 
ducible. We set w{a a) = {) for notational simplicity. 
The escape rate \{a) is then defined by 



A(a) = w{a — ^ a'). 



(Gl) 



a'eM 



Let P''^{a,t) be a probability of a{t) — a, where the 
superscript w represents the transition rate dependence 
of the probability P™. Then, P'^{a, t) satisfies the master 
equation 



d_ 

di 



a' 

(G2) 

For later convenience, we define a matrix by 

A{a,a') =w{a' ^ a) - \{a)5a,a' ■ (G3) 
The stationary distribution P^{a) satisfies 

5]^(a,a')P,T(a')=0. (G4) 
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We denote by (g)^ the expectation value of a function 
g{a) with respect to the stationary distribution P^{a). 

We fix a time interval [0,r]. For a given history 
{a{t))l^Q, where a{t) £ A4, a. sequence of transition times 
{ti,t2, ■ ■ ■ ,tn~i) is determined, and a{t) is expressed as 
a(t) = ai for ti-i < t < U with to — and t„ = t. 
Here, we consider a quantity a{ai, ai-^-i) defined for the 
jump event — > a^+i ai t = U {1 < i < n — 1). Its 
time-averaged value in the history is written as 



^ 71 — 1 

A{t) = - a{ai,ai+i). 



(G5) 



We then study the scaled cumulant generating function 
of A(t)^ which is defined by 



G(K) = lim — log ( e 

r—^OQ T 



ThA{r) 



(G6) 



where ( ) represents the sample average over realizations 
of history with an initial distribution. 

It is known that Gili) is equal to the largest eigenvalue 
of the matrix An given by 



Ahip! , a) = w{a — > a')e 



l\ ha{a,,a') 



Ka)5a,a'. (G7) 



See Ref. [SJ] for the derivation. Let 0max(a) be the left 
eigenvector for the largest eigenvalue of Ah- That is, 



^(j)m!,y,{a')Ah{a' ,a) = G(h)(j)^^^{a). 



(G8) 



Since (/)max(a) > (See Appendix [BJ , we may define a 
potential function T4nax(a) as 

Vmax(a) = -21og0max(a). (G9) 

Dividing the both sides of (jG8[) by 0max(ci), we obtain 
G{h) 

= ^ \u]{a a')e-V'max(a')/2+y„ax(a)/2+ftQ(a,a') 
a' 

— w{a -> a') . 

(GIO) 

Here, for any potential function V{a), we define a modi- 
fied transition rate u\[a — > a') as 



uX,{a a') 

= w{a a')e-^(''')/2+^(°)/2+''"("^"') 
The corresponding escape rate A)f becomes 

a' 

and (jGlOp is rewritten as 

G{h)^X];--{a)-Xia). 



(Gil) 



(G12) 



This expression might be known, although we do not find 
its explicit presentation in previous studies. 

Now, our main claim, which is far from trivial, is the 
variational expression 



G{h)=max{X]:-X)2\ (G14) 



and the relation 
dG{h) 



dh 



Y,P:t' {a)u]r{a^a')a{a,a') (G15) 



with 



Vopt = argmax (A^ - A)^ 



(G16) 



The derivation is the following. First, we consider the 
variation V ^ V + 5V m. ^A^ — A)^^** , which is expressed 



as 



5P:J (a) [XX (a) - A(a)] + J] (a)^^ (a). 

a a 

(G17) 



From (jGlip and (|G12p . we rewrite the second term of 
(IGTtI) as 

a a' 

(G18) 

Since this becomes zero from the condition of 
the stationary distribution, the variational equation 



leads to 



K = XUa)-Xia), 



(G19) 



where if is a constant corresponding to a Lagrange mul- 
tiplier. Here, (|G19I) is equivalent to 



)/2 



(G20) 



By comparing this with (IG8I) . we find that V^opt = Knax 
as the result of the Perron-Frobenius theory. 

Next, we consider the modified system with the tran- 
sition rate it^°'''(a — > a'). By multiplying the both sides 

of (|G13p by Pgj'' (a) and summing over a, we obtain 



G(;i) = (A^'-A)^ 



(G21) 



where we have used Vopt = Knax- We thus conclude 
(|G14[) . Furthermore, by differentiating (IGISP with re- 

(G13) spect to h, multiplying the both sides by P^^.'^ (a), and 
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summing over a, we derive 
dG{h) 



dh 



a, a' 



dh 
(G22) 



The second term of (jG22[) turns out to be zero from the 



condition of the stationary distribution. The final ex- 
pression is (jGlSp . 



At the end of this Appendix, we point out that the 
expression (jG14[) is similar to the formula of the large 
deviation function of the empirical measure in Markov 
jump processes !553. This similarity corresponds to the 
relation between ([7^ and (15^ in Langevin equations. 
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